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Abstract We briefly review sources of cosmic rays, their composition and 
spectra as well as their propagation in the galactic and extragalactic mag- 
netic fields, both regular and fluctuating. A special attention is paid to the 
recent results of the X-ray and gamma-ray observations that shed light on 
the origin of the galactic cosmic rays and the challenging results of Pierre 
Auger Observatory on the ultra high energy cosmic rays. The perspectives 
(-h i of both high energy astrophysics and cosmic-ray astronomy to identify the 

O || sources of ultra high energy cosmic rays, the mechanisms of particle accelera- 

tion, to measure the intergalactic radiation fields and to reveal the structure 
£~ ' of magnetic fields of very different scales are outlined. 
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1 Introduction 



To a first approximation, the all-particle spectrum of cosmic rays (CR) - 
ionized atomic nuclei with energies extending from the MeV to more than 
10 20 eV, can be described by a power law over 11 decades on particle en- 
ergy, so that the dependence of cosmic ray intensity on particle energy is 
close to E~ 2 - 7 above 10 GeV/nucleon see Figured] Spectra of different CR 
species are different however. Closer examination reveals some structures in 
the galactic cosmic ray spectrum that includes more flat low energy part at 
E < 10 GeV/nucleon (mainly due to the modulation of the Galactic cosmic 
ray intensity by the solar wind flow), the knee at about 4 x 10 15 eV, the 
possible second knee at ~ 4 x 10 17 eV, the ankle at ~ 4 x 10 18 eV where the 
spectrum flattens again, and finally the suppression of flux above ~5x 10 19 
eV. The latte r can be interpreted as the so-called Greis en-Zatsepin-Kuzmin 
(GZK) cutoff (jGre isen 1966l lZatsepm and Kuzmin1ll966t) caused by the inter- 
action of the ultra-high-energy cosmic rays (UHECRs) with the cosmological 
microwave background (CMB). Therefore, the suppression of the UHECR 
flux, ta ken together with the detection of an anisotropy in their arrival di- 
rection ( Abraham et al.l [20071 ) and the absence of any anisotropy associated 



with galactic structure is often taken as evidence of the extragalactic origin 
of these particles. However, it should be noted that the reported cutoff at 
10 20 eV could also be related to the limited capability of even extreme cosmic 
accelerators to boost the energy of individual particles up to 10 20 cV. In any 
case, it is widely believed that the hig hest energy cosmic rays above 10 19 eV 
are of extragalactic origin (see however iBlasi et al]l200dl"Calvez et al]l20ich . 

On the other hand, the energy region below 10 15 eV is dominated by 
particles produced by Galactic sources. A strong argument for this is pro- 
vided by gamma-ray observations of the Galactic disk, and the compari- 
son with the corresponding observations in the Magellanic clouds. Further- 
more, it is likely that even low-intensity magnetic fields in the intergalac- 
tic space prevent low-energy CRs from reaching us from distant galaxies 
on cosmological timescales. Therefore, a transition from galactic to extra- 
galactic cosmic rays is expected to occur somewhere between 10 17 eV and 
10 19 eV, while it is tempting to associate the flattening of the spectrum 
at the so-called ankle in the energy spectrum as marking the transition re- 
gion. However, di fferent interpretations have been proposed, with a transitio n 
at lower energy ( Berezinskv et al.l 120041 . iLemoinel 120051 lAloisio et alJl2007h . 



and th is important part of high-energy CR phenomenology is still much de- 
bated (|Allard et al.ll2007llBerezinskvll2009l ). It is interesting to note that this 



question is strongly related to that of the source composition of extragalac- 
tic CRs (lAUard et aUl2005h . and of course, the energy scale of the transition 



depends on the capacity of the Galaxy to confine cosmic rays up to high 
energies, and thus on the intensity and structure of the Galactic magnetic 
field. 

The spectra of both Galactic and extragalactic cosmic rays are shaped by 
two basic processes - the acceleration in the sources and the subsequent prop- 
agation in cosmic magnetic fields and radiation fields. We concentrate here 
only on the CR propagation and acceleration while the extensive review of 
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Fig. 1 All particle spectrum of cosmic rays (after S. Swordy and from Abbasi ct al. 
(20081) ~ insert). 



the galactic and extragalacti c magnetic fie l ds is presented in the other chap- 
ters_of the volume (see also lFerrierdl200ll iBeckl [20081 iKulsrud and Zweibell 



I2008L and the references therein) 



2 Potential sources 

If one uses the cosmic ray energy requirements and the nonthermal radiation 
as a guideline, then the most powerful accelerators of rclativistic particles in 
the Galaxy should be supernovae and supernova remnants, pulsars, neutron 
stars in close binary systems, winds of young massive stars, and possibly 
gamma-ray bursts or hypcrnovae. The total power L cr needed to maintain 
the observed energy density of cosmic rays is estimated as 10 41 erg/s. For 
the acceleration by a supernovae, this estimate implies the release of energy 
in the form of cosmic rays of approximately 1.5 x 10 50 erg per supernova if the 
supernova rate in the Galaxy is 1 every 50 years. This value comes to about 
10-20% of the kinetic energy of the ejecta, which is compatible with the 
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expectations of the theory of diffusive shock acceleration. This assumes the 
acceleration of cosmic rays by the outward propagating shock, which results 
from the supernova explosion and propagates in the interstellar medium or 
in the wind of the progenitor star. 

The rotational energy of a young pulsar with period P that remains af- 
ter the supernova explosion is estimated to be 2 x 10 50 (10 ms/P) 2 erg. It is 
an additional energy reservoir for particle acceleration and in particular it 
could be the source of very high energy electron-positron pairs. The produc- 
tion of relati vistic pos i trons by pulsars may explain their prese nce in c osmic 
rays, see e.g. lProfumd(l2008h . A general dimensional analysis bv lHillasI (|l984l 
120051 ) provided a useful tool to estimate the maximal CR energies that can 
be achieved by electromagnetic acceleration of charged particles in a source 
of given size and magnetic field magnitude. It follows that the highest en- 
ergy CR particles of energy > 10 20 eV can occur either in the compact 
sources with very large magnetic fields or in the most extended cosmological 
objects, though the particle energy losses may affect the maximal CR energy 
estimations. 



2.1 Supernova remnants as CR accelerators and gamma-ray astronomy 

The observations of non-thermal radio emission of supernova remnants have 
established the pr esence of ef f ective acceleration of cosm ic ray electrons in 
these objects, e.g. lLozinskaval ( 1992t ). iDrurv et al.1 ( 200ll ). Shell- type super- 



nova remnants exhibit a broad range of spectral indices, centered roughly on 
a = —0.5. This implies a power-law distribution of electrons with an average 
index 7 = 2a + 1 w 2.0, but it varies in a broad range from about 1.4 to 
2.6. The an alysis of the synch rotron emission in the young supernova rem- 
nant Cas A ( Jones et al.ll2003h showed the presence of electrons with energies 



up to 200 GeV with a magnetic field strength of about 500 fiG. The inter- 
pretation of nonthermal radio emission from external galaxies proved that 
supernova remnants are the sites of acceleration of relativistic electrons with 
an efficiency similar to the on e needed to provi de the observed intensity of 
Galactic cosmic-ray electrons ([Duric et al.lll995l ). 

In the case of young SNRs, the spectra of nonthermal synchrotron radia- 
tion continues into the X-ray region. In particular, bright nonthermal X-ray 
structures have been observed from SN 1006, RX J1713. 7-3946, RCW 86, 
and Cas A. The only viable explanation of this emission is the synchrotron 
radiation of ultrarelati vistic electrons accelerated t o energies as high as 100 
TeV (for a review, see IVinkl [200l iRevnoldsl 120081) . Inverse Compton scat- 
tering of background photons by these electrons, together with gamma rays 
produced via ir° decays are two major mechanisms of gamma-ray production 
related to the forward and reverse shocks. 

It is generally believed that gamma-ray astronomy has a key role to play 
in solving the problem of origin of galactic cosmic rays. The basic concept is 
simple and concerns both the acceleration and propagation aspects of cosmic 
rays. Namely, while the localized gamma-ray sources exhibit the sites of pro- 
duction/acceleration of cosmic rays, the angular and spectral distributions 
of the diffuse galactic gamma-ray emission provide information about the 
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character of propagation of relativistic particles in galactic magnetic fields. 
The crucial outcome of current efforts involving both the ground-based and 
space-borne gamma-ray detectors is hoped to be a decisive test of the popular 
hypothesis that supernova remnants (SNRs) are responsible for the bulk of 
the observed cosmic ray flux. Three factors, (i) high efficiency of transforma- 
tion of the available kinetic energy of bulk motion into nonthermal particles, 
(ii) the predicted, within the DSA paradigm, hard energy spectra of protons 
extending to multi-TeV energies and (iii) relatively high density of the ambi- 
ent gas - make the young supernova remnants viable sources of gamma-rays 
resulting from production and prompt decay of secondary 7r°-mcsons. Thus, 
the most straightforward test of acceleration of proton and nuclei in SNRs 
can be performed via search for 7r°-decay gamma-rays - either directly from 
shells of young SNRs ( Drurv et al.lll994f) or from ne arby clouds interacting 



with an expanding SNR shell ( Aharonian et al.lll994|) 



One may argue that the TeV gamma-ray domain is the best energy band 
to explore this possibility - from the point of view of the superior performance 
of the detection technique and because of the key information about particle 
acceleration carried by TeV gamma-rays. On the other hand, TeV gamma- 
rays are expected only from young SNRs, when the shock speed can be as 
high as 3000 km/s - a key condition for acceleration of protons to energies 
100 TeV and beyond. With the age, typically t > 2000 yr, the shock speed 
and the magnetic field in the expanding shell are decreased to the point 
when the particle acceleration significantly slows down. Moreover, the highest 
energy particles accelerated at the early epochs, effectively leave the remnant. 
This results in dramatic reduction of gamma-ray emissivity at TeV energies. 
Meanwhile the production of GeV gamma-rays continues. Thus we expect 
MeV and GeV gamma-rays not only from young, but, to a large extent, 
also from mid-age SNRs. In this regard we expect significantly more GeV 
emitting SNRs compared to the TeV emitting SNRs. Unfortunately, due to 
relatively low gamma-ray fluxes, as well as very high gamma-ray background 
in the crowded regions of the galactic disk, the detection of GeV gamma- 
rays from SNRs is quite difficult. In the case of mid-age SNRs we expect also 
" delayed" gamma-ray emission produced by cosmic rays which have left the 
source and interact with nearby massive gas molecular clouds. The spectrum 
of cosmic rays, and, consequently, gamma-rays produced in p-p interactions 
depends on several factors, including the spectrum of protons that escaped 
the supernova remnant, the character of diffusion, the age of the source and 
the distance between the source and cloud. In this regard one should expect 
gamma-ray emission with GeV-to-TeV flux ratio which can vary significantly 
from a source to source. 

Over the last several years TeV gamma-ray emission has ben reported 
from several SNRs. Presently six young shell type SNRs - Cas A, RX J1713.7- 
3946, RX J0852-4622 (Vela Jr), RCW 86, SN 1006, Tycho - are identified as 
TeV gamma-ray emitters. Several others, in particular IC 433, W28 and 
CTB 37B are, most likely, SNR/molecular-cloud interacting systems. Also, 
several galactic TeV gamma-ray sources spatially coincide with the so-called 
composite supernova remnants, objects with characteristic features of both 
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stand ard shell-type SNRs and Pulsar Wind Nebulae (see e.g. lAharonian et al.1 
12008ft . 

In Fig|3] we shown the TeV and X-ray images of the strongest and best 
studied gamma-ray emitting SNR - RX J1713. 7-3946. 




Fig. 2 The gamma-ray image of RX J1713. 7-3946 obtained with the H.E.S.S. tele- 
scope array (|Aharonian et aljr2007bf ). In the left-hand side, the overlaid light-gray 
contours illustrate the significance of the different gamma-ray features at 8, 18, and 
24a levels. In the right-hand side, the black-lines show the X-ray contours. 

The age of this SNR RX J1713. 7-3946 at a distance of about 1 kpc is 
estimated between 1 to 3 thousand years. Therefore, it can be formally teated 
as a representative of young galactic SNRs like SN 1006, Vela Junior, Tycho, 
Kepler, Cas A. However, RX J1713. 7-3946 seems to be a unique object with 
quite unusual characteristics. This circumstance makes the identification of 
the gamma-ray mechanism, and hence the nature of the parent particles a 
nontrivial problem. Meanwhile the question of the origin of gamma-rays has 
a fundamental implication related to the role of SNRs to the production of 
galactic cosmic rays. 

Remarkably, both the hadronic (gamma-rays produced in pp interactions) 
and leptonic models (gamma-ray produced via inverse Compton scattering) 
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provide, with more or less success, but yet acceptable fits to the measured 
gamma-ray spectra and explain also the spatia l distributions of gamma-rays 
and X-ravs IZirakashvili and Aharonianl ( 2010|) . So the test of these models 



relies essentially on multiwavelength data. 

One of the most puzzling features of RX J1713. 7-3946 is the lack of 
thermal X-ray emission which is a serious argument against the hadronic 
models of of this source. The tight upper limit on the thermal X-ray flux of 
RX J1713. 7-3946 is explained by the supernova explosion inside a wind-blown 
bubble with a very low gas density n gas <C 1 cm~ 3 . This allows t o suppress the 



free- free component of thermal emission but not the X-ray lines ([Ellison et al.1 
120101 IZirakashvili and Aharonianl2010l ). thus an additional assumption about 
the peculiar composition with reduced co ntent of heavy elements is require d 



to avoid the conflict with X-ray data ( Zirakashvili and Aharonianl l2010l) . 



The hadronic models demand quite large magnetic field of order of 100 fiG 
or larger. This, for the given flux of synchrotron radiation, requires only 
10 46 erg in electrons. On the other hand, the total energy in protons 
should be 10 50 erg or significantly larger if the background gas density n <C 
1 cm~ 3 . Thus, in the hadronic scenario the electron-to-proton ratio is close to 
K ep = 10~ 4 . This is a two orders of magnitude smaller than the "standard" 
electron/proton ratio observed in local cosmic rays. 

A major challenge for the leptonic models is the demand for a small 
magnetic field in the shocked shell of the remnant. Within the framework of 
a one-zone model in which the electron acceleration and gamma-ray emission 
take place in the same region, the magnetic field cannot significantly exceed 
10 /iG, otherwise it leads, for the given X-ray flux, to significant suppression 
of IC gamma-rays. 

A possible solution to the problems related to the pure (oversimplified) 
one-zone leptonic and hadronic models could be effective acceleration of elec- 
trons and protons in both forward and reverse shocks with an additional as- 
sumption of existence of dense condensations embedded in very low density 
shell. In this scenario the overall VHE gamma-ray emission is dominated by 
the IC scattering of electrons, but in addition to the leptonic component one 
should expect non-negligible contribution from hadronic gamma-rays pro- 
duced at interactions of protons with dense gas condensations. Moreover, in 
this scenario both the low energy, E < 100 GcV and highest, E > 10 GcV 
gamma-rays should arrive from specific regions occupied by dense gas con- 
densation (see FigEJ) ■ The "composite" scenario allows a quite relaxed pa- 
rameter space compared to the "pure" hadronic model and leptonic models, 
and at the same time, provides better fits to the broad-band spectral energy 
distribution of the source. 

Independent of the ability of different models to describe the spectral 
and morphological features of gamma-ray emission of RX J1713. 7-3946, it is 
obvious that we deal with a source that effectively accelerates electrons and 
protons to energy of 100 TcV and beyond. On the other hand, this object 
has many unique features, and in this regard it could be misleading if we 
treat this source as a representative of the whole population of young SNRs. 
Therefore, the studies of other TeV emitting SNRs is crucial for definite 
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Fig. 3 Broad-band emission of RX J1713. 7-3946 for the "composite" scenario of 
gam ma-rays with non-modified forward and reverse shocks and dense clouds in the 
shell Zirakashvili and Aharonian (2010). 



conclusions concerning the role of SNRs in production of galactic cosmic 
rays. 

The realization of different possible scenarios of particle acceleration and 
radiation in general, and the uncertainties in relevant model parameters like 
the strength of the magnetic field and the density of ambient gas, prevent 
us from unambiguous conclusions concerning the fraction of the hadronic 
component in the flux of observed gamma-rays, and consequently do not al- 
low model-independent estimates of thew total energy released in accelerated 
protons and nuclei. The best hope in this regard is related to the detection 
of high energy neutrinos produced in interactions of accelerated protons and 
nuclei with the ambient plasma. Although the fluxes of neutrinos arc quite 
similar to the gamma-ray fluxes, because of limited sensitivity of the neutrino 
detectors, this extremely important channel of information unfortunately re- 
mains unexplored. Even for the largest neutrino telescopes like the recently 
completed IceCube and the planned several km-cube underwater detector in 
the Mediterranean Sea, the high energy neutrino signal from RXJ 1713.7- 
3946 is expected to be marginal even under the extreme assumption that 
the flux of dete cted TeV gamma-rays is fully consists of photons of hadronic 
origin (see e.g. IVissani et al.l [20TTI iMorlino et al.l [20091 and the references 
therein). 



2.2 OB-associations, starburst, supcrbubblcs 



Collective stellar winds and supernovae in clusters and associations of mas- 
sive sta rs in galaxies can be pot ential ly efficient cosmic ray accelerators 
(see e.g. Cesarskv and Montmerle 19821 Bvkov and Fleishmanlll992L Axforc 
1994LlHigdon et al.lll998llBvkovll200lLlBvkov and ToDtvginll200lLlBinns et al 



20071) . 

The young rich open cluster We sterlund 2 is spatially c oincident with 
TeV photon source HESS J1023-575 (|Aharonian et al.ll2007al) . The reported 
source extension excludes a single star origin of the observed VHE emission 
and argues in favor of the production of cosmic rays by collective stellar 
winds from the ensemble of at least a dozen of hot massive O-type and a 
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few WR-type stars in the stellar cluster Westerlund 2 (of an estimated age 
of about a few million years). 

Recently, TcV gamma-ray emission ha s be reported from two star-burst 
galaxies, NGC 253 and M82, by H.E.S.S. (lAcero et al.ll2009h and VERITAS 
( Acciari et alJl2009h . respectively. Starburst galaxies exhibit in their central 
regions a highly increased rate of supcrnovae that are likely to be efficient 
producers of cosmic rays. As they are extended, with radii up to ~ 100 pc, 
the superbubblcs (SB) containing dozens of supernova remnants and winds 
of massive stars arc difficult to identify in the Galaxy, and are thus currently 
best studied in the Magellanic Clouds. Nonthermal X-rays were detected by 
iBamba et all (|2004 from the shells of the superbubble 30 Dor C in the LMC. 
The X-ray observations revealed the SB morphology as a nearly circular 
shell with a radius of 40 pc, which is bright on the northern and western 
sides. The non-thermal X-ray shell traces the outer boundary of the radio 
shell. These features of thin thermal and non-thermal X-rays are similar to 
those of SN 1006, a prototype of a synchrotron X-ray shell, but the non- 
thermal component of 30 Dor C is about 10 times brighter than that of SN 
1006. The source is also much older than that of SN 1006, so the particle 
acceleration time in this superbubble may be longer than those in normal 
shell-like supernova remnants. 

Nonlinear modeling of CR acceleration in superbubbles bv iBvkovl ( 200ll ) 
predicted a high efficiency of the kinetic energy conversion to the CRs that 
can be ~ 20% and even higher with a significant temporal evolution of particle 
spectra inside the superbubble. Magnetic field inside the superbubble can be 
amplified up to about 30 fiG providing the maximal energies of CR protons to 
be about 10 s GeV and higher for CR nuclei that dominate the accelerated CR 
composition at highest energies in the model fsee lBvkov and To ptvein 200 lh . 
To model the CR composition in superbubbles a number of components 
i ncluding the 22 Ne-rich Wolf-Rayet and dust material mu st be accounted for 
( Binns et al ] l2007l iMever et al.l(l997l lEllison et al.lll997h . 



2.3 CR abundance constraints on potential CR sources 

The correction of the cosmic ray composition observed at the Earth for the 
effects of nuclear fragmentation in the interstellar medium makes it possible 
to determine the initial elemental and isotopic composition of accelerated 
particles, to clarify the process of cosmic ray acceleration and the nature of 
cosmic ray sources. 

The relative abundance of chemical elements in cosmic ray sources is in 
general similar to the solar and local galactic abundances but with some 
interesting deviations. The popular scenarios which explain the cosmic ray 
source composition include the acceleration of grains together w ith relatively 



less abundant volatile ion s by shocks in the interstellar medium (jMever et al.1 
119971 [Ell ison et al"1ll997j ). the acceleration of freshly formed materia l, partic- 
ularly grains in young supernova remnants ( Lingenfelter et aT1ll998[ ). and the 



acceleration in ho t superbubbles with multiple supernova remnants (see e.g. 
iBinns et al.ll2007h . 
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The isotopic composition of cosmic rays is now measured for all stable 
isotopes for elements from H to Ni. The isotopic composition of cosmic ray 
sources is strikingly sim ilar to the composition found in the solar system 



nysrmi 

(|Wiedenbeck et allboOlh . However, a well-established anomaly in the iso- 



topic composition of galactic cosmic rays is the excess of 22 Ne. The ratio 
22 Ne/ 20 Ne is enhanced by a factor of 4 compared with the solar reference 
value (jDuvernois et al.lll996aD . It can be explained only by the special con- 
ditions of nucleosynthesis. The enhancement of neutron rich isotopes would 
be expected in the highly evolved very massive stars in their Wolf-Rayet 
stage when their surfaces cont ain large excesses of th e products of core he- 
lium burning, including 22 Ne ( Casse and PaulHl982f ). An increased cosmic 



ray ratio C/O by a factor of 2 is also in favor of Wolf-Rayet stars. 

Using the data from Cos mic Ray Isotope Spectrometer aboard the ACE 
spacecraft iBinns et al.l (]2007t ) measured the isotopic abundances of neon and 
a number of other species in the galactic cosmic rays. The authors found 
that the three largest deviations of galactic cosmic ray isotope ratios from 
solar-system ratios, namely 12 C/ 16 , 22 Ne/ 20 Ne, and 58 Fe/ 56 Fe, are consis- 
tent with a model of CR source consisting of about 20% of Wolf-Rayet wind 
material mixed with about 80% material of the solar-system composition. 
Therefore [Binns et al.l (|2007t) concluded that OB associations within super- 



bubbles are the likely source of at least a substantial fraction of GCRs. 

The analysis of the 59 Co/ 56 Fe abundance ratio also provides important 
i nformation. Ind eed, the relatively high, close to the solar value, of this ratio 
(jMewaldtlfl999h testifies that the major part of the originally synthesized 



59 Ni, which is stable when it is ionized (i.e. when it is accelerated among 
the cosmic rays), has decayed by the K-capture of an orbital electron into 
59 Co before the acceleration started. Thus, the delay between the synthesis 
of this material and its acceleration has to be larger than 10 5 yr severely 
constraining the models with the acceleration of freshly ejected matter in 
SNRs. 

The production of cosmic rays with power L cr = 10 erg/s per normal 
galaxy can provide a cosmic-ray energy density of WMg = L cr N g Tn = 3 x 
10~ 17 erg cm -3 , where N g = 2 x 10 -2 Mpc -3 is the number density of 
galaxies and Tjj = 1.4 x 10 10 yr is the Hubble time. This is almost 5 orders 
of magnitude smaller than the energy density of cosmic rays observed in the 
solar system. The contribution of active galaxies can reduce this discrepancy 
by only a factor of 10. All this is in agreement with the picture described 
in the Introduction, where the main part of the cosmic rays we observe at 
the Earth is produced by galactic sources (e.g. supernova remnants) and 
slowly diffuse in galactic magnetic fields out of the Galaxy to the intcrgalactic 
space, where the cosmic ray density is much smaller. The efficiency of the 
confinement of energetic particles is decreasing with energy, and because 
of that the difference between the cosmic ray densities inside and outside 
the Galaxy is also decreasing. The densities become equal at some energy 
between 10 17 and 10 19 eV and the extragalactic sources dominate at higher 
energies. The list of potential extragalactic sources includes in particular 
the active galactic nuclei, the progenitors of gamma -ray bursts, the accretion 
shocks of the size of galactic clusters and larger fsee lTorres and Anchordooul 
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|2004 ISigJ I2009L for a review). These objects are probably able to provide 
the required power q^ g = 3 x 10 36 erg s _1 Mpc -3 in cosmic rays with 
energies > 10 19 5 eV. This estimate of qu a is based on the cosmic ray flux 
determ ined using the Auger Observatory (|Pierre Auger Collaboration et all 
l2010aTl and it takes into account the energy loss time due to the GZK effect. It 
is notable that the entire popul ation of supernovae i n the Universe also meet 
the required power condition (lAublin and Parizotl 120061) . but as far as we 
know the highest energy particles can not be accelerated in such sources. The 
origin of the highest energy CRs in extragalactic sources remains a mystery, 
whose solution will require enhanced detections capabilities (because of the 
very low fluxes involved, see below), but will have a strong impact on the 
whole field of high-energy astrophysics. 

3 Transport in the Galaxy 

The interaction of relativistic charged particles with galactic magnetic fields 
explains the high isotropy and relatively large confinement time of the cosmic 
rays in the Galaxy. It is accepted that the diffusive approximation gives 
an adequate description of the cosmic ray propagation in the interstellar 
medium at energies up to ~ 10 17 eV. The diffusion model forms a basis 
for interpretation of the cosmic ray data and the related radio-astronomical, 
X-ray and gamma-ray observations. 

3.1 Empirical model of cosmic ray propagation 

The procedure to model the cosmic ray diffusion in the Galaxy can be sum- 
marized as follows. One must first specify the cosmic ray sources, define 
the shape of the cosmic ray halo and the conditions at its boundaries (it 
is usually assumed that energetic particles exit freely into the intergalac- 
tic space where the cosmic ray density is negligible). The basic diffusion- 
convection equations for different cosmic ray species describe the diffusion 
of relativistic charged particles and, if needed, their convective transport 
in the models including a galactic wind. The equations should incorporate 
all possible energy loss and gain processes in the interstellar medium, nu- 
clear fragmentation, and the radioactive decay of unstable nuclei. One can 
then calculate the distribution functions of protons, electrons and the dif- 
ferent types of nuclei. The transport coefficients of cosmic rays (diffusion 
coefficient and convection velocity), the properties of the cosmic ray sources 
(total power, energy spectra of the different components, elemental and iso- 
topic composition), and the size of the confinement region of cosmic rays in 
the Galaxy can be inferred by fitting all the available data on cosmic rays. 
Hundreds of stable and radioactive isotopes are included in the calculations 
of nuclear fragmentation and the transformation of the energetic nuclei as 
they interact with the interstellar gas. The most advanced code developed 
for the numerical modeling of cosmic ray propagation in the Galaxy is the 
GALPROP code which uses a Crank-Nicholson im plicit second-order scheme 
([Strong and Moskalenkolll998l IStrong et al.l l2007h . It incorporates as much 
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realistic astrophysical input as possible, together with the latest theoreti- 
cal development, and numerically solves the transport diffusion-convection 
equations for all cosmic-ray species. 

One of the key channels of information about cosmic-ray propagation is 
the abundance of the secondary energ etic nuclei Li, Be, B , Sc, V, Ti, 2 H, 
3 He and others produced as a result of the spallation of heavier primary nu- 
clei interacting with the interstellar gas. The observed flux ratios between 
secondary and primary nuclei, for example the boro n to carbon ratio , show 
a maximum at about E = 1 GeV /nucleon (e.g., ( Engelmann et al.l ll990L 
lYanasak~ ct al. 200ll lAhn et al.f2008f ). The mass distribution of secondary nu- 
clei can be understood in terms of their propagation through about 10 g/cm 2 
of interstellar material. This quantity, referred to as the escape length, can be 
approximately expressed through the parameters of the plain diffusion model 
with a flat cosmic ray halo of total thickness 2H pa 8 kpc as X ps V[i g H/2D, 
where v is the particle velocity, /i 9 « 2.4 mg/cm 2 is the surface gas density of 
galactic disk, and D is the cosmic ray diffusion coefficient. This allows an es- 
timate of the diffusion coefficient: D ~ 3 x 10 28 cm 2 /s at E = 1 GeV/nucleon. 

The primary-to-secondary ratios decrease at higher energies (at least up 
to a few hundred GeV/nucleon where the measurements exist) as well as at 
lower energies. The understanding of this behavior is provided by the kinetic 
theory of particle interaction with the interstellar magnetic fields. 



3.2 Kinetic theory of cosmic ray diffusion in galactic magnetic fields 

The theory of energetic particle transport in the galactic magnetic fields is 
constructed in much the same way as in the better studied case of particle 
transport in the interplanetary magnetic fields. The detailed tre atment of 
cosmic ray diffusion in magnetic fields can be fo und in m onographs (|TqptvgiiJ 
ll985l . lBerezinskii et alJll990USchlickeiserl2002l ) (see also lPtuskin et alfeOOaT 
On the "microscopic level" the spatial and momentum diffusion of cosmic 
rays results from the particle scattering on random MHD waves and disconti- 
nuities. The effective "collision integral" for energetic charged particles mov- 
ing in small-amplitude random fields SB <C B can be taken from the standard 
quasi- linear theory of plasma turbulence, see e.g. iKennel and Engelmann! 
(119661) . The wave-particle interaction is of resonant character and depends 
on the particle gyroradius, which is equal to r g = pc/ZeB = \Q 12 Rqv cm in 
a typical interstellar magnetic field B = 3 fiG (here Rqv is the particle mag- 
netic rigidity in GV and [i is the cosine of particle pitch angle). The particle 
is predominantly scattered by the irregularities of magnetic field SB that 
have the projection of the wave vector on the magnetic field direction equal 
to feii = ±s/ (r g /i) . The integers s = 0,1,2... correspond to the cyclotron 
resonances of different orders. The efficiency of particle scattering depends 
on the polarization of the waves and on their distribution in k-space. The 
first-order resonance s = 1 is the most important for the isotropic and also for 
the one-dimensional distribution of random MHD waves along the average 
magnetic field, k || B. In some cases - for calculation of scattering at small /i 
and for calculation of perpendicular diffusion - the broadening of resonances 
and magnetic mirroring effects should be taken into account. Locally, the 
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cosmic ray diffusion is anisotropic and goes along the local magnetic field. 
The isotropization of diffusion is accounted for by the presence of strong 
large-scale (~ 100 pc) fluctuations of the galactic magnetic field. The prob- 
lem of calculation of the average diffusion tensor is not trivial even in the case 
of relatively weak large scale random fields, since the field is almost static 
and the strictly one-dimensional diffusion along the mag netic field lines doe s 
not lead to non-ze ro diffusion perpendicular to B, see ICasse et al.l ( 20021) . 
IWebb et all (120091) . 

The following simple equations for the diffusion coefficient are useful for 
estimates: 



vrl B 2 
3(l-a)kl- a 5B~l- 



D = vr g B 2 / [127rk les W(k les )} = — ^7757 , (1) 



here fc rcs = l/r g is the resonant wave number, and W{k) is the spectral 
energy density of waves normalized as J dkW(k) = 5B 2 /4w. The random 
field at the resonance scale is assumed to be weak, 5B rcs <C B. The principle 
wave number of the turbulence is k^ and the amplitude of random field at this 
scale is SB^. The second equation in (1) is valid for particles with r g < fc^ 1 
under the assumption of a power law spectrum of turbulence W(k) oc 1/ fc 2_a , 
k > fc L . The diffusion coefficient has then a power law scaling on momentum 
D oc v{p/Z) a . 

Eq. (JTJ) can be used to describe the turbulence which consists of random 
isotropic distribution of Alfven and fast magnetosonic waves and it gives 
correct order of magnitude estimates for the wave distribution concentrated 
around the direction of average magnetic field. It should be pointed out 
however that the isotropization of the diffusion tensor does not occur in the 
case of a pure parallel wave propagation (k || B). Another special case is 2D 
turbulence with perpendicular propagation of waves (k 1 B). In this case, 
the scattering occurs only for magnetosonic waves through the resonance 
s = which leads to a very large diffusion coefficient, about a factor {v/V a ) 2 
larger than given by eq. ([I]). 

The diffusion in momentum is approximately described by the equation: 

D pp =p 2 V 2 /(9D). (2) 

Here equal intensities of waves moving along the magnetic field in opposite 
directions are assumed. The momentum diffusion is vanished {D pp = 0) if all 
waves move in one direction along the magnetic field with the same phase 
velocity. The flow of waves carry particles along the field in this case and the 
convection with velocity determined by the phase velocity of waves appears 
in the transport equation for cosmic rays, see Skilling (1975a), Berczinsky et 
al. (1990) for detail. Clearly the cosmic ray convection also arises as a result 
of the large scale motion of the interstellar gas with frozen magnetic field. 

The random component of interstellar magnetic field with an extended 
spectrum of inhomogencities can give the resonant particle scattering and 
spatial diffusion of cosmic rays. Information on the interstellar turbulence 
spectrum have been obtained from radio scintillation and refraction obser- 
vations (sensitive to fluctuations of thermal electron density) , measurements 
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of the differential Faraday rotation angles from distant sources (mainly pro- 
duced by fluctuations in the interstellar magnetic field), and the observations 
of random motions in the interstellar gas. These data are consistent with the 
assumption that a singl e close-to-Kolmogorov s pectrum extends from scales 
10 8 to 3 x 10 20 cm, see lArmstrong et al.1 (|l995t) and references therein. The 
Kolmogorov spectrum is of the form W(k) oc fc~ 5 / 3 . Other types of spectra 
frequently used to describe the interstellar turbulence are W(k) oc k~ 2 for 
the sh ock-dominated turbulence, see e.g. the model by Bvkov and Toptyginl 
(11987), an d spect rum W(k) oc fc~ 3 / 2 suggested bv llroshnikovl ( l*964h and 
iKraichnanl ( 1965t ) in their phenomenological theory of MHD turbulence. 
Comprehensive reviews of MHD tu rbulence with appl ic ation to the interstel- 
l ar co nditio ns have been given by IZhou et al.l (|2004D .IeI megreen and Scalol 
(|2004h . and lScalo and Elmegreenf (j2004h . 

The estimate based on the empirical value of the diffusion coefficient for 
GeV particles (see section 13. 1[) gives the level of turbulence at the princi- 
pal scale &l = 1CP 21 cm -1 of the order 8B tot /B ~ 0.2 for an Iroshnikov- 
Kraichnan spectrum W(k) oc fc -3 / 2 (a = 1/2) and SB to t/B ~ 1 for a 
Kolmogorov-type spectrum W(k) oc fc~ 5 / 3 (a = 1/3). At the same time, 
the data on Faraday rotation angles favor the Kolmogorov spectrum with 
SBtot/B ~ 1 and k^ = 10~ 21 cm" . The cascades of Alfven waves (with the 
scaling fc~ 5 / 3 ) and the fast magnetosonic waves (fc~ 3 / 2 ) are independent and 
may coexist in the iGoldreich and Sri dharl (119951) model of MHD turbulence. 
The amplitude of Alfven wave cascade may dominate at the principle scale. 

There are two physical explanations of the observed secondary-to-primary 
nuclei ratios in the diffusion model with a static cosmic ray halo. 

The first explanation ( Simon et al.lll98(j ISeo and Ptuskinlll994[) refers to 
the diffusion model with distributed reacceleration of cosmic rays by the in- 
terstellar MHD turbulence which scatters particles and provides their spatial 
diffusion. The Kolmogorov-type spectrum is assumed that leads to the rising 
with rigidity of the diffusion coefficient D oc v(p/Z) 1 ^ 3 . For typical value of 
Alfven velocity V a ~ 30 km s _1 , the reacceleration is not essential for nuclei 
with energies E > 40 GeV/nucleon and the abundance of secondary nuclei is 
a decreasing function of rigidity. The impact of reacceleration on spectra of 
primary and secondary nuclei becomes stronger at smaller energies so that the 
characteristic time of distributed acceleration in the Galaxy becomes equal 
to the time of diffusion from the Galaxy at about 1 GeV/nucleon. In con- 
sequence the pronounced peak in the secondary-to-primary ratio arises, see 
Figure ((H). The asymptotic behavior of the escape length is A" oc v(p/Z)~ 1 ^ 3 
at E > 40 GeV/nucleon. 

The second explanation ( Ptuskin et al.l2006l) assumes that the Iroshnikov- 
Kraichnan spectrum describes the interstellar MHD turbulence. It is charac- 
terized in particular by a relatively slow nonlinear cascade of wave from small 
to large wave numbers. The resonant wave-particle interaction results in the 
significant wave damping on cosmic rays and termination of the cascade at 
k ~ 10 -12 cm -1 (in contrast to the case of Kolmogorov cascade which is fast 
and not noticeably affected by cosmic rays). The amplitude of short waves 
is suppressed and the low energy particles rapidly exit the Galaxy without 
producing many secondaries. It explains the peaks in secondary/primary nu- 
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Fig. 4 B/C ratio as calculated in plain diffusion model with unphysical break of 
the diffusion coefficient (a), reacceleration model (a), and diffusive reacceleration 
with damping model (c). Lower curve - local interstellar spectra, upper - modu- 
lated (with modulation parameter <& = 450 MV). Data below 200 MeV/n ucleon: 
ACE rtPavis et alj l200Ct ) . Ulysses (|Duvernois, Simpson, anTT havcr 19963), Voy- 
ager (ILuka siak 199% ; high energy data: H EAO-3 (lEngelmann et al.H1990l), C REAM 
(|Ahn et al.N2008i ) for other references see IStephens and Streitmatterl 1)19981 ). 



clei ratios at about 1 GeV/nucleon, see Figure ffl. The asymptotic behavior 
of the scape length in this case is X oc u(p/Z) -1 / 2 at E > 3 GeV/nucleon. 

It is clear from Figure (U) that the experimental data on the abundance of 
secondary nuclei are well described by both physical models. Also illustrated 
in Figure (j?]) are the results of calculations in the plain diffusion model with 
the empirical diffusion coefficient of the form D oc /3~ 2 and f3~ 2 (p/Z) - 6 
below/above 3 GV chosen to fit the B/C data (the "unphysical" model). 

It must be emphasized that the theoretical description of MHD turbulence 
is a complicated and not completely solved problem even in the case of small- 
amplitude random fields. Since the mid 1990s, there has been a renewed inter- 
est in understanding of magnctohydrodynamic turbulence as it applies to the 
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interstellar magnetic field an d density fluctuations ( Goldreich and Sridharl 



ll995lll997llGaltier et al.ll200(Th . lGoldreich and Sridharl (| 1995h exploited anisotropy 



in MHD turbulence and obtained Kolmogorov-like spectrum for the energy 
density of Alfvcn waves. The main part of the energy density in this turbu- 
lence is concentrated perpendicular to the local magnetic field wave vectors 

1 /2 

k± ~ k, while the parallel wave numbers are small: fc|| ~ [kW(k)/ (Sq/47t)] k±. 

2/3 

The cascade is anisotropic with energy confined within the cone fc|| oc kj_ . 

This turbulence does not significantly scatter cosmic rays. The distribution of 

slow magnetosonic waves passively follows that of Alfven waves ( Lithwick and Goldreich! 



l200ll ). They are not damped on cosmic rays (because of the property k± ^> 
fc||) and are probably responsible for the observed interstellar electron density 
fluctuations. The fast magnetosonic waves with the Iroshnikov - Kraichnan 
spectrum W(k) oc fc -3 / 2 may have an independent nonlinear cascade which 
is isotropic and can efficiently scatter cosmic rays. These concl usions were 
supported by numerical simulations by ICho and Lazarianl ( 2002h . This con 



cept of the MHD turbulence favors the second of the discussed above sce- 
narios where cosmic rays are scattered by fast magnetosonic waves with the 
Iroshnikov-Kraichnan spectrum and suppress the turbulence at large wave 
numbers k > 10 -12 cm' 1 . 

The knowledge of the diffusion coefficient is absolutely essential for un- 
derstanding the nature of the spectrum of galactic cosmic rays that is deter- 
mined by the processes of acceleration in the sources (supernova remnants) 
and propagation in galactic magnetic fields. Two specific asymptotic power 
laws of the diffusion coefficient D oc (p/Z) 1 / 3 and D oc (p/Z) 1 / 2 at very high 
energies together with the observed spectrum approximated by the power 
law J oc E~ 2 - 7 imply the cosmic ray source spectra close to q oc E 22 and 
q oc E 2A respectively. 

The contemporary modelling of cosmic ray production by supernova rem- 
nants was achieved bv lPtuskin et al.l (|2010f ). The spectra of high-energy pro- 
tons and nuclei accelerated by supernova remnant shocks were calculated 
taking into account magnetic field amplification and Alfvenic drift both up- 
stream and downstream of the shock for different types of supernova rem- 
nants during their evolution. Four different types of SNRs (la, Ib/c, IIP and 
lib) with corresponding burst rates were included in the calculations. The 
action of cosmic ray pressure on the shock structure was taken into account 
in the calculations. It was found that the maximum energy of accelerated 
particles may reach 5 • 10 18 eV for Fe ions in Type lib SNRs. The steady 
state spectrum of cosmic rays produced by SNRs in the Galaxy was calcu- 
lated with the deduced source spectra of protons and different kind of nuclei 
up to Iron. Th e escape leng t h with a high energy asymptotic X oc (p/Z)~°- 54: 
determined by I Jones et al.l ( 200ll ) by the accurate fit to B/C data was used 



in the calculations. The derived energy spectrum of cosmic rays including the 
knee structure around 3 x 10 15 eV is in good agreement with the spectrum 
measured at the Earth from low energies to about 5 x 10 18 eV. This result 
is strongly in favor of the cosmic ray galactic diffusion on the MHD turbu- 
lence with spectrum close to k~ 3 l 2 . The Kolmogorov type spectrum k~ 5/>4 is 
unlikely. 
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The diffusion approximation cannot be used when the diffusion mean free 
path I = 3D I v exceeds the size of cosmic ray halo H w 4 kpc that occurs at 
E > 2 x 10 16 Z eV if D oc (p/Z) 1 / 2 . Calculations of particle trajectories in 
galactic magnetic field arc needed to study cosmic ray propagation at these 
ultra high energies. At weaker dependence, D oc (p/Z) 1 / 3 , the diffusion ap- 
proximation breaks at higher energies. The main limitation in this case is the 
inefficiency of particle scattering at r g k^ > 1 when the diffusion coefficient 
along the average magnetic field rapidly rises with energy as D ~ w^l, see 
e.g lToptvginl (jl985t ). This regime is realized at E > Z x 10 17 eV. The direct 
numerical mode lling of particle moti on in the regular and random galactic 
magnetic fields ( Horandel c t al. 2007|) confirmed the presented order of mag- 
nitude estimates. The most essential feature of cosmic ray transport at ultra- 
high energies above about 10 17 eV is the prevalence of the drift motion across 
predo minantly azimuthal average galactic magnetic field ( Zirakashvili ct al. 
Il998l ). It leads to the inverse dependence of the exit time from the Galaxy 
on particle momentum T oc (j>jZ)~ x . 



3.3 Instabilities and plasma effects 

Cosmic rays not always can be tr eated as test particles moving in given 



Cosmic rays not always can be treate d as test particles moving m given 
regular and random mag netic fields (|Ginzburdll9MlWentzellll974LICesarskvl 



I1980L lKulsrudll2005l ). The energy density of cosmic rays estimated aS UUcr — 
1 — 2 cV cm -3 is approximately equal to the energy density of the magnetic 
field and to the energy density of turbulent motions of the interstellar gas. 
The presence of relativistic charged particles in the interstellar medium leads 
to collective (plasma) effects during cosmic ray propagation. 

A notable example is the cosmic ray streaming instability which de- 
velops when the bulk velocity of cosmic rays exceeds the Alfven velocity 
V a . The growth rate of MHD waves amplified by relativistic charged parti- 
cles with number density n cr (E) oc i?~ 7+1 at the resonance wave number 
feres = r~ 1 (E) is 

where f2 p is the gyrofrequency of thermal protons, S CT is the amplitude of 
cosmic ray anisotropy. Even for small anisotropy, <5 cr ~ 10 -3 , the instabil- 
ity for galactic cosmic rays with energies ~ 100 GeV develops in about 10 5 
years, i.e. rather rapidly for the galactic timescale. The development of the 
instability leads to isotropisat ion of th e angu l ar distribution of partic l es and 
turbulence amplification, see IZweibell ( 2003|) . iFarmer and Goldreichl ( 2004 ) 



and references cited therein. Plasma effects make the overall picture of cosmic 
ray diffusion in the Galaxy more complicated than it was discussed in Sec- 
tion 3. In principle, the cosmic ray diffusion coefficient should be calculated 
sclfconsistcntly with the account taken for the generation of turbulence by 
streaming cosmic rays. The examples of such approa ch and the corre spond- 
ing transport e quatio ns can be found in the papers bv lSkillingl (1975a.b) and 



iPtuskin et al.1 (|2008h . The necessity of considering the variety of dissipation 
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processes of linear and nonlinear wave dissipation in the interstellar medium 
adds complexity to the investigation and makes the results strongly depen- 
dent on the interstellar gas parameters: the density, the state of ionization, 
the temperature. The wave dissipation makes the effect of streaming insta- 
bility not efficient at high energies since the cosmic ray density is going down 
with particle energy and the growth rate Eq. [3] is decreasing correspondingly. 
One of the exceptions is the model of galactic wind dri ven by cosmic rays in a 
rotating galaxy developed bv lZirakashvili etahl ((l996) (it is briefly described 
at the end of the present Section) . Moving predominantly along the very long 
spiral magnetic held lines, the cosmic rays have enough time to amplify the 
resonant waves even at high energies of the order 10 15 eV. 

The effect of the streaming instability is more significant in the vicinity of 
the sources. The amplification of magnetic field by the cosmic-ray streaming 
instability at the shock in supernova remnant is an integral part of cosmic 
ray acceleration (lBelj|l978l iBlandford and Eichlerll 19871 iMalkov and Drur 
20011 lBellll2004 IZirakashvili and Ptuskinll2008l iMarcowith and CasseH2oT 



Bvkov et al. 201 ll ). The cosmic ravs produce turbulence which selfconsis 



tently determines the cosmic ray diffusion coefficient in the shock vicinity. 

Cosmic rays may induce, in addition to kinetic effects, significant hydro- 
dynamic effects in the Galaxy. Accounting for cosmic ray pressure is princi- 
pally important for the formation of a halo filled with gas, magnetic field, and 
rclativistic particles. The equilibrium distribution of the interstellar medium 
above the galactic plane in the gravitational field of stars is subject to the 
Parker instability ( Parker! fl966l ). Cosmic rays play a significant role in the 



development of this instability. The instability gives rise to large-scale tur- 
bulence and helps sustain an almost cquipartition energy distribution among 
cosmic rays, magnetic fields, and turbulent gas motions. The characteristic 
time for instability development is ~ 10 7 year s in the gaseous galactic disk, 
and ~ 10 8 years in the gas halo. iParkerl ( 19921 ) showed that magnetic arches 



and loops appearing above the galactic disk due to the action of cosmic rays 
are necessary for au dynamo to operate, which is primary mechanism of 
magnetic field generation in the Galaxy. A numerical model of the magncto- 
hydrodynamical d ynamo, driven by co smic rays in the interstellar medium 
was developed bv lHanasz et al.l (|2004) . The cosmic rays in the model were 
accelerated in randomly occurring supernova remnants. The cosmic ray prop- 
agation was accounted with t he diffusion-advectio n equation supplementing 
the MHD equations (see also lHanasz et aill2006l ). The other essential ele- 
ments of the model are: vertical gravity of the disk, differential rotation and 
resistivity leading to reconnection of magnetic field lines. They obtain ampli- 
fication of the large-scale magnetic field on a timescale of galactic rotation. 
The model i ncludes the id eas of a fast, cosmic ray driven, galactic dynamo 
proposed bv lParke"rl ( 1992ft . The authors find that both the resistivity and the 



introduced SN activity enhance the efficiency of the cosmic-ray dynamo. The 
times cale of magnetic fie ld amplification in this model was as short as 140 
Myr (jHanasz et al.|[2006l ). The models assumed some CR diffusion prescrip- 
tions. A linear study of the effect of the anisotropic CR diffusion with very 
different b ut finite parallel {k\\) and transverse diffusion coefficients was per- 
formed bv lRvu et al.1 (|2003ft . It has been shown that the finiteness of parallel 
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Fig. 5 The structure of the galactic wind flow driven by cosmic rays l|Ptuskin et al.l 
119971 ). The large-scale magnetic field lines, the flow streamlines and the trajectories 
of cosmic ray bulk motion from the Galaxy are confined in the same flux tubes. The 
cylindrical flux tube with the cross section A(s) originating at the Sun location is 
shown (here s is the coordinate along the flux tube). Diffusion mainly determines 
the cosmic-ray transport in the diffusion zone (the analog of the boundary layer) , 
whereas convection mainly makes it in the advection zone. The boundary between 
two zones s*(p) is moving up with energy of cosmic-ray particles. 



diffusion slows down the development of the Parker instability. However, with 
a realistic value of km = 3 x 10 28 cm 2 s~ x in the ISM, the maximum growth 
rate is smaller by only a couple percent than that for km — > oo, and the range 
of unstable wavenumbers remains the same. The inclusion of perpendicular 
diffusion with Kn = 0.02k^ docs not change the growth rate noticeably. It 
should be noted however, that the turbulent advection due to magnetic fluc- 
tuations of scales larger than the CR mean free path would make the diffusion 
coefficient nearly isotropic on scales interesting for the Parker instability (see 
iBvkov and Toptvginlll993l). A non-linear development of Parker instability 
was simulated by iKuwabara et all ()2004[ ). They stated that the growth rate 
of the instability is larger in the models with smaller diffusion coefficient only 
in the early linear stage. On the later stages the growth rate becomes smaller 
when compared to that of the large diffusion coefficient model. The growth 
of instability is impeded by the CR pressure gradient force interfering with 
the falling motion of the matter in the small- km model, while the magnetic 
loop can grow up to larger scales in the models with large K||. 
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It is possible that the gas in the galactic halo is not in static equilib- 
rium but is involved in large-scale convectivc motion - the galactic wind. 
The data o n galactic soft X-ray emission suggest that a wind exists in our 
Gala xy, see Everett et al] (12007). It can be supported by cosmic ray pres- 



sure (llpavich 19751 Breitschwerdt et alJll99ll ). A model was constructed by 
IZirakashvili et al.1 (|l996f ) where cosmic ravs. after leaving the sources (super- 
nova remnants), determine the wind outflow in the rotating Galaxy with a 
frozen magnetic field, see Figure ([3]). Here, the streaming instability of cos- 
mic rays exiting the Galaxy along the spiral magnetic field lines leads to the 
MHD turbulence generation, that self-consistcntly determines the diffusion- 
convection transport of relativistic particles. The level of turbulence is regu- 
lated by the nonlinear Landau damping on thermal ions. The outflow velocity 
is ~ 30 km s _1 at a distance ~ 20 kpc, and the speeds up to a velocity of 
~ 400 km s _1 several hundred kpc away. The external pressure of the inter- 
galactic gas produces a termination shock at a distance of ~ 300 kpc. In this 
model, the diffusion coefficient of cosmic ray is not given independently and 
is self-consistently calculated, being d ependent on the pow er of sources and 
the spectrum of accelerated particles ( Ptuskin et al.llT997l ). Remarkably, the 
obtained transport coefficients and other model parameters are consis tent 
with the empirical diffusion-convection model, e.g. lBloemen et ahl ( 19931 ). for 
cosmic ray propagation in the model with galactic wind. It should be em- 
phasized that the existence of a galactic wi nd, either the r mal o r driven by 
cosmic rays, remains hypothetical. However. lEverett et al.l ( 2008f ) found that 
the observed galactic soft X-ray emission can be better explained by a wind 
than by the static gas models and the cosmic ray pressure is essential to drive 
the wind. 



4 The highest energy cosmic rays 

The " Hillas diagram" discussed a bove predicts th at the compact sources like 
gamma ray bursts (G RB) (see e.g. lMeszaroslj2006l. for a review) are potential 
sources of UH ECRs (|Vietrilll995l IWaxmanlll995l iMilgrom and Usov|[l995l 
IWaxmar] 12004 ) . as well as the cos mological sources like activ e galactic nu- 
clei (AGN), radi o galaxies (see e.g. Longairll2010L lBlasill2006l ). and clusters 



of galaxies (e.g. iNorman et al.1 Il995h . The observed widespread warm-hot 
intergalactic gas was likely heated by large scale cosmological shocks. The 
shocks in the Large Scale Struc ture formation sc enario are driven by gravita- 
tionally accelerated flows (e.g. iKang et"al 1 I2005D . Non-thermal emission ob- 
served from clusters of galaxies indicate the presence of acceler ated particles 
in scales larger the galactic halo sizes fe.g. iFerrari et al1l2008|) . Large-scale 
accretion shocks in the cluster of galaxies and more generally cosmological 
shocks in Large Scale Structure can a ccelerate particles to energies above 
10 19 eV. Namely, H orman et al ] (119951) argued that cosmological shocks can 



be good sites for UHECR acceleration if there is an intergalactic field of 
order 10~ 9 G and even more important if microgauss regime fields can be 
self-generated in shocks accelerating particles. The maximal energies cosmic 
rays accelerated in the large scale shocks like the cluster accretion shocks or 
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hot spots in radio galaxies can reach ~ 10 19 eV reg ime only if indeed mi- 
crogauss fields can be p roduced (iNorman et al. f995). Curren t cosmological 
shock simulations (e.g. iBvkoy et al.l 20081 Dolaget al.ll2008[ ) and the non- 
thermal emission models (e.g. Kushnir et alj 20091 ) are in favor of the reality 
of microgauss regime magnetic fields in clusters of galaxies. On the other 
hand the losses caused by interactions of the UHECRs with photons of the 
Cosmic Microwave Background Radiation could limi t the maximal energ ies 
of protons accelerated in clusters of galaxies (see e.g. lVannoni et aT1l2009l) . 



4.1 Ultra- high energy cosmic rays as probes of magnetic fields 

At the highest energies, the gyroradius of the cosmic rays becomes extremely 
large, especially in the intergalactic space where the magnetic fields are ex- 
pected to be much lower than in the Galaxy. Their mode of propagation thus 
changes from being diffusive to being asymptotically rectilinear, so at extreme 
energies it is realistic to think of the possibility of a cosmic ray astronomy, 
with the sources of ultra-high energy cosmic rays (UHECRs) being identi- 
fiable directly. For example a proton of energy 6 10 19 eV will be deflected 
by 1° to 5° in the galactic magnetic field depending upon the direction and 
length of the trajectory. 

On the one hand, being able to actually see the UHECR sources would 
be a remarkable breakthrough: i) it would (at least partially) solve the mys- 
tery of their origin, ii) it would put strong constraints on the acceleration 
mechanism and its efficiency, teaching us a lot about the physics of parti- 
cle acceleration in the universe in its most extreme manifestations, and iii) 
it would offer an important piece of information, complementary to that 
gathered from multi- wavelengths observations of photons, to guide us in the 
astrophysical modeling of the sources. This would indeed enhance the power 
of current multi-messenger analysis, which are mostly relying on upper limits. 

On the other hand, if we could be sure to know the source of (at least 
a set of) UHECRs, we could use these particles as messengers, not only to 
gather new information about the sources, but also to probe the intervening 
magnetic fields, by studying the deflections they suffered along their journey 
to the Earth. 

In the following, we review the current status of UHECR observations, 
and discuss some of their implications. We also discuss why a new generation 
of detectors can be expected to lead to the above-mentioned breakthrough, 
and what can already be learned from the current data, including from the 
fact that its statistics is still insufficient. But to understand the data, we first 
need to understand the transport of UHECRs, and see how it can modify 
their energy spectrum, their composition and their arrival directions. 



4.2 UHECR transport: energy, mass and direction 

The key ingredient of UHECR transport in the extragalactic space is the 
so-called GZK effect, which refers to the interactions of the hi gh-energy cos- 
mic rays with the background radiation fields in the universe (|Greisenlll966l 
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IZatsepin and Kuzminl [l966h . In the case of high-energy protons, the domi- 
nant interactions are with the photons of the cosmological microwave back- 
ground (CMB), producing electron-positron pairs and/or pions as soon as 
the proton energy is high enough for the CMB photons to exceed the reac- 
tion threshold, when boosted in the proton rest frame. These interactions 
lead to strong energy losses, especially above ~ 6 10 19 eV, when photopion 
production becomes important, which signific antly limits the distance UHE 
protons can propagate in the universe (see e.g. iBhattachariee and Sigl|[2000l 
for a review). This translates into an effective horizon, from beyond which 
extragalactic sources cannot contribute to the observed UHE protons, and 
this in turn implies an expected strong suppression of the flux, referred to as 
the GZK cutoff. At lower energy, energy losses due to photopair-production 
on the CMB photons also modify the spectrum of UHE protons, produc- 
ing a dip around 3 10* 8 eV. where the ankle is observed in the overall CR 
spectrum ( Blumenthallll970llBerezinskii and Grigor'evalll988l ). Another 



sim- 



ple interpretation of the ankle, however, is that it marks the transition from 
galactic to extragalactic cosmic rays. 

As already noted in the original work of lGreisenl ( 19661) and lZatsepin and Kuzminl 



(|l966f ). a similar effect holds for heavier nuclei, which can be photodissociated 
by CMB and/or infrared photons as they propagate through the intergalactic 
space. The propagation of ultra-high - energy nuclei in the intergalactic space 
was studied in detail bv lPuget et al. (fl976l). and recently revisited with up- 



dated and additional cross-sections ( Khan et al.ll2005l lAllard et al.ll2005l ). A 



key result is that photodissociation is very efficient at ultra-high-energies, 
and lead to: i) a modification of the composition of UHECRs as they prop- 
agate from their sources to the Earth, as well as ii) energy losses associated 
with the loss of nucleons (although pair production also plays a role at lower 
energy, to a first approximation the nuclei essentially keep a constant Lorentz 
factor). 

The energy threshold for the photodissociation of UHE nuclei depends on 
the energy of the background photons, as well as on the mass of the nuclei. 
At the threshold, the dominant process is the excitation of the giant dipolc 
resonance, involving photons around 10 MeV in the nucleus rest frame. Pho- 
tons in the tail of the CMB distribution reach this energy for nuclei with a 
Lorentz factor of the order of -TJ>10 9 , and energies E~4x I^m p c 2 , where A 
is the mass number of the nucleus. Thus, heavier nuclei can survive photodis- 
sociation up to higher energies. By a striking coincidence of Nature, although 
the processes involved are very different, and the cross section thresholds and 
sizes are also different, the energy loss length for UHE protons and for Fe 
nuclei are essentially the same in the intergalactic space, so both species have 
a very similar horizon structure. Therefore, their elemental GZK cutoff leads 
to the same spectral shape, and indeed the energy distribution of UHECRs 
can be understood just as well in terms of pure proton sources or of pure Fe 
sources (or a combination of both, provided the source spectral index and/or 
the evolution of the source luminosity as a function of redshift is adjusted 
correspondingly). In principle, intermediate mass nuclei may also be accel- 
erated in the same UHE sources. However, at a given energy, their Lorentz 
factor is smaller than the Lorentz factor of Fe nuclei, so they see higher en- 
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ergy photons and get dissociated more rapidly. As a consequence, the very 
highest energy cosmic rays are likely to be mostly protons and/or Fe nuclei 
(heavier nuclei would also survive, but a priori they are much less abundant 
in the sources). 

The studies on the propagation of UHE nuclei by lAllard et al.l (|2005l 

120071 ) also drew a renewed attention to the question of the galactic/extragalactic 
transition, whose energy scale turns out to depend on the assumed UHECR 
source composition. Schematically, a pure proton composition (at the sources) 
would allow one to fit the high energy spectrum down to below 10 18 eV, as- 
suming a steep source spectrum (in E~ a , with a ~ 2.6 or 2.7, depending 
on the assumptions on the cosmological source evolution), and would then 
imply a galactic/extragalactic transition at a few times 10 17 eV. Conversely, 
a mixed source composition (e.g. similar to that of low energy cosmic-rays in 
the Galaxy) would allow one to fit the energy spectrum with a source spectral 
index a ~ 2.3, and imply that the extragala ctic component overcomes the 
galactic one at the ankle ( Allard et alJl2005l) . We note in passing that this 
question is also related to the question of the galactic magnetic fields, which 
must be able to confine the cosmic rays up to the transition energy 

Any information about the UHECR composition would thus have impor- 
tant implications for the general phenomenology of cosmic rays, including 
low energy ones. Unfortunately, a corollary of the above remarks is that the 
overall UHE spectrum cannot be used, on its own, to constrain the source 
composition of the UHECRs. A direct measurement of the mass of UHE par- 
ticles is thus very important, and although this is experimentally difficult, 
noticeable progress has been made recently (see below). However, it is fair 
to say that the experimental situation is not clear yet, and since the mea- 
surements also rely on hadronic physics at energies well beyond the energy 
accessible through direct experiments, it is not clear how conclusively this 
kind of measurements can resolve the issue in the near future. 

In this context, two other aspects of UHECR transport can play a role, 
which are worth mentioning here since they also involve magnetic fields. First, 
the interactions associated with the above-mentioned GZK effect (s) produce 
secondary particles, including the so-called cosmogenic neutrinos and the 
electromagnetic cascades initiated by secondary gamma-rays (from pion de- 
cay) or electrons and positrons. These cascades involve the usual mecha- 
nisms of synchrotron radiation, pair production and i nverse Compton interac- 
tions (see e.g., in the context of UHECR propagation. lBhattachariee and Sigl 
l2000t ) . These secondary particles can in principle be detected (or their fluxes 
constrained). In the case of secondary neutrinos, the expected flux and en- 
ergy distribution strongly depends on assumptions about the phenomenol- 
ogy of the UHECR sourc es (for a recent update, including sce narios with 
a mixed composition, see lAllard et al.l 120061 iKotera et al.l 20101 and refer- 
ences therein). It should also be noted that secondary particles produced in 
the acceleration site may be detected or constrained, and thus provide ad- 
ditional, multi-messenger information about the UHECR sources and their 
environment, as well as about the acceleration mechanism itself. Finally, ex- 
tragalactic magnetic fields can also be constrained by the obser vation of TeV 
gamma-ray halos around potential acceleration sites. Recently Tavlor et al.l 
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(2011) discussed the constraints on the strength of extragalactic magnetic 
field derived from the simultaneous detections of TeV emitting blazars in 
GeV photons by the Fermi telescope. The measured GeV flux was found to 
be lower than that calculated in the cascade model under the assumption of 
zero magnetic field. Assuming that the reason for the suppression of the cas- 
cade component is the extended nature of the cascade emission, the authors 
concluded that the extragalactic magnetic field above 10 -15 G of correlation 
length of 1 Mpc is consistent with the data. 

The second aspect, directly associated with magnetic fields, is the de- 
flection of UHE particles between their sources and the Earth. If the source 
of each individual cosmic ray was known, the angular distance between the 
source and the arrival direction would give a very precious measurement of 
the deflection integrated along the path of the cosmic ray. Since the deflection, 
in a ballistic or semi ballistic regime, is proportional to the particle charge, 
an analysis of the deflection patterns would provide information about the 
UHECR charges and/or about the intervening magnetic fields. Although the 
identification of sources apparently remains out of reach of the current de- 
tectors, we shall mention below what the current observational situation can 
nevertheless teach us, and why this situation is expected to change with a 
ten times larger statistics at the highest energies. It should also be mentioned 
that additional hints about the composition of UHECRs might come from 
the study of their anisotropics, comparing patterns in the arrival directions at 
different energies and/or angular scale. The basic idea is that protons of en- 
ergy E have the same gyroradius as heavier nuclei with energy Z x E. Now, 
in a regime where energy losses can be neglected in a first approximation 
(say, below the GZK energy scale), the angular transport of the different 
UHE particles only depend on their gyroradii. Therefore Fe nuclei from a 
given source follow exactly the same path as protons at an energy 26 times 
lower. This has recent l y bee n turned into a very interesting argument by 
iLemoine and Waxmanl (|2009t ). relating anisotropy studies at different ener- 
gies to composition assertions. Extensive studies with larger statistics should 
provide important clues about this crucial issue in the future. Likewise, at the 
highest energies, if the number of sources is low and a few hot spots appear 
in the UHECR sky, halos of different angular sizes around the centroid of the 
source should appear, corresponding to different nuclear species (essentially 
protons and Fe or sub-Fe nuclei, if both are present at the source). 

More generally, it is clear that the study of the UHECR angular deflec- 
tions have a lot to teach us, both about the cosmic rays and about the inter- 
vening magnetic fields. The main aspects of this transport can be summa- 
rized as follows. We may distinguish deflections occurring: i) at or around the 
sources, ii) in the intergalactic space, and iii) in the Milky Way. Even though 
the sources remain unknown, it is not unlikely that they are embedded in 
environments where large scale magnetic fields exist, e.g. in galactic clusters. 
These magnetic fields may be able to deflect even the highest energy CRs 
significantly, or even confine them for some time before they escape from the 
source environment. This can have an effect of the energy spectrum reaching 
out in the intergalactic space, if the confinement time becomes comparable 
to the energy loss time. It may also affect the composition of the UHECR 
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leaving the sources, since, at a given energy, heavier nuclei are more easily 
confined. However, this would require effective lengths of several tens of Mpc 
at ~ 10 20 eV. While such a confinement would increase the angular size of the 
source on the sky, turning it, say, from a localized source to the size of a whole 
cluster of galaxies for instance, this should not affect our ability to identify 
clusters of events from individual sources on the sky, as long as the number 
of sources is limited (as expected at the highest energies) and their typical 
angular separation is larger than the apparent size of the confinement halo. 
Finally, it is worth noting that a specificity of UHE cosmic rays, compared 
to the low-energy, galactic ones, is that the constancy of the fluxes cannot 
be assumed anymore, since the dynamical timescales of the sources can be 
shorter than the confinement timescales around them and/or the differen- 
tial propagation timescales to the Earth. This may result in time dependent 
effects (on time scales of typically thousands to tens of thousands of years) 
as to which sources contribute to the observed flux at different energies, as 
well as in different mass ranges, depending on the elapsed time since their 
activity as UHECR sources. 

Once they leave their sources, UHECRs may also be deflected by mag- 
netic fields in the intergalactic voids, if their strength is large enough. This, 
however, is highly uncertain. Among the interesting effects that may be as- 
sociated with such magnetic fields, the existence of a magnetic horizon has 
received particular attention, as it was realized that it could act as an anti- 
GZK effect, i.e. essentially a high-pass filter for energetic particles, prevent- 
ing lower energy cosmic rays from reaching us from distant sources, if they 
wander around in the extragalactic space (due to magnetic diffusion) for a 
time longer than the energy loss timescale, or even the a ge of the universe 
(|Parizotli200ilLemoinel[2005llAloisio and Berezinsky||2005l) . This can have an 
impact on the matching between the galactic and extragalactic components, 
and also on the cosmic-ray composition in the corresponding energy range: 
at a given energy, heavier nuclei diffuse for a longer time than light ones, 
and thus extragalactic Fe nuclei should not be able to reach us from distant 
galaxies at a energy as low as protons do. It should be noted, however, that 
the actual outcome of such propagation effects on the U HECR spectrum and 
composition depend on the granularity of the sources (|Delignv et alj|2004 
lAloisio and Berezinskvl |2005|) . If the sources were uniformly distributed in 
the universe, no effect at all would be expected, so no definite prediction can 
be made at our current stage of ignorance about the sources, the source den- 
sity and of course the intergalactic magnetic fields. More recent works have 
focused on the pr opagation of UHECRs in s trongly inhomogeneous extra- 
galactic magnetic ( Kotera and Lemoinell2008c] lal) . This is indeed a much more 
plausible assumption, under which UHECRs appear to propagate essentially 
in straight lines between more localized interactions with large scale magne- 
tized structures, acting as scattering centers. As a consequence, potential hot 
spots in the UHECR sky may reveal the position of the most nearby scatter- 
ing cen ters rather than that of the most nearby sources ( Kotera and Lemoind 
l2008al) . 

Different aspects of prop agation of cosmic rays in IGMF r ecent l y have 

been comprehensively studied lDolag et all (l2005l ). lGlobus et all (|2008l ). lKotera and Lemoind 
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Fig. 6 Energy distributions of protons observed within different angles for the 
source at distances r — 100 Mpc (left panel) and r = 300 Mpc (right panel). The 
energy spectrum of protons is assumed in the form of power-law with exponential 
cutoff, with a = 2 and E = 3 x 10 20 eV; the rate L = 10 44 erg/s. 



(|2008bl ). The conclusions of these studies are generally different, basically be- 
cause of different assumptions and approaches in the modeling of the IGMF. 
The propagations of protons and acc ompanying them gamm a-rays and neu- 
trinos recently have been studied bv lAharonian et "ail ( 2010h in small angle 
approximation limit. In Fig|5] we show the energy distributions of protons 
within different angles arriving from a point source located at distances 100 
and 300 Mpc. It is assumed that the source is injects protons into the inter- 
galactic medium with magnetic field InG and correlation length A = 1 Mpc. 
The upper dashed lines in Fig(5] correspond to the case when protons propa- 
gate in empty space; flux is determined by the geometrical factor 1/r 2 . The 
solid line present the case when the deflections in the magnetic field are ig- 
nored. Comparison of these two curves demonstrates strong dependence of 
the energy distribution of protons on the solid angle within which they are 
dctct-ttHl: the flux of protons at highest energies is concentrated along the 
direction to the source. This is the result of selective deflection of of protons 
depending on their energy. In addition to this effect there are additional two 
spectral features - a bump and a sagging at lower energies. The bump pre- 
ceding the cutoff is due to strong growth of energy losses at the threshold 
of photomcson production that makes particles to be accumulated in this 
energy region; the sagging is a consequence of the energy losses due to the 
electron-positron pair production. However, as it is seen in FigjSJ the maxi- 
mum due to the deflection is more distinct than the "photoproduction energy 
loss bump" . 

Finally, at the other end of the propagation track, the galactic magnetic 
fields also deflect the UHECRs, with deflections of the order of a degree times 
the charge of the particle, or possibly much more if an extended magnetized 
wind surrounds our Galaxy. However, the exact deflection patterns (both size 
and direction) depend on the structure of the magnetic field, which remains 
highly uncertain. The regular component should provide a global shift on the 
sky of the image of the actual UHECR sources, while the turbulent compo- 
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nent essentially broadens the apparent cosmic-ray source around that image. 
Finally, the deflections are expected to be inversely proportional to energy, 
which should lead to characteristic structures in the arrival direction/energy 
space. Chains of UHECR events, ordered in decreasing energy away from 
what would be a common source, have been sought in the various data sets, 
with no conclusive result so far. Obviously, while such structures would be 
extremely informative about the source positions as well as about the mag- 
netic field lines responsible, they can be detected only if the sources do not 
overlap (i.e., at the highest energies, but then, unfortunately, with a limited 
lever arm in energy), or with much larger statistics, to allow the detection 
of these patterns on top of a background mixing UHECRs from overlapping 
sources, with a large enough significance. 

In conclusion, the study of UHECR transport in the energy space, mass 
space and angular space, makes it clear that the road to a better understand- 
ing of the sources and phenomenology of UHECRs passes through the joint 
study of the three fundamental obscrvables (energy spectrum, composition 
and arrival directions). It should also be stressed that the GZK cutoff, while 
implying a limitation of the UHECR energies in the universe (independently 
of the sources) and strongly reducing the expected flux, is also the best chance 
we have to answer the central questions about UHECRs. As the current data 
show (see below), no obvious source and/or pattern has appeared in the sky 
yet. This implies either that the sources are numerous, and we have not seen 
multiplcts (i.e. several events coming from the same astrophysical source), 
or not strong enough for angular associations in the sky to be significant, 
or that the overall deflections of UHECRs are large, and the sources largely 
overlap. The obvious way out is by reducing the number of sources as well 
as the deflections. This is precisely what the GZK allows (and actually im- 
poses). By reducing the horizon scale as the energy increases, the GZK effect 
guarantees that fewer and fewer sources contribute at the highest energies. In 
addition, the number of sources capable of accelerating cosmic rays up to a 
given energy is likely to be a decreasing function of that energy. In addition, 
the deflections are inversely proportional to energy, so by concentrating on 
the UHECR sky at 10 20 eV or so, we are bound to face a situation where a 
(very) limited number of sources contribute to the observed flux, with limited 
deflections and well separated "centroids" of the corresponding hot spots in 
the UHECR sky. The price to pay for such a situation, of course, is that the 
corresponding cosmic ray fluxes are very low, so further progress is condi- 
tioned by the resolution on an important observational challenge. The target 
has been set at the level of a total exposure of the sky of several 10 5 up t o 
10 6 km sr yr (see the recent white papers on UHECRs: lOlinto et alJl2009ft . 
i.e. more than an order of magnitude larger than what is currently available. 

4.3 UHECR Observations 

The highest-energy particles are so rare that they are detectable only by 
means of the giant cascades or extensive air showers they create in the at- 
mosphere. Details of how these extensive air showers are observed and of 
how the parameters of importance are measured can be found in reviews 



28 



(e.g. iNaeano a nd Watson 2000). The results that are most relevant to the 
topic of this chapter have been obtained with the Pierre Auger Observa - 
tory (|Abraham et all |2004 IT) awson and for the Auger Collaboration! |2008[ ) 
in Argentina and with the HiRes instrument ( Bover et al. 12002^ operated in 
the Dugway Desert, Utah, USA. The HiRes device uses telescopes compris- 
ing mirrors and photomultipliers to detect the fluorescence radiation emitted 
from the excitation of atmospheric nitrogen by the shower particles passing 
through the atmosphere. The Auger Observatory has similar telescopes but in 
addition has an array of water- Cherenkov detectors deployed over 3000 km 2 
to measure the particles in the showers at ground level. The fluorescence 
telescopes are used to obtain a measure of how the cascade grows and decays 
in the atmosphere while the water-Chercnkov detectors are used to obtain 
the pattern of the distribution of shower particles on the ground from which 
a shower size, closely proportional to the primary energy, is obtained. With 
both devices the direction of the incoming cosmic ray is measured with an 
accuracy of ~ 1° using GPS-based timing. 

A 21 st century innovation has been to build instruments combining the 
fluorescence and the surface detector techniques. This combination has al- 
lowed what is known as the hybrid method to be developed in which the 
time of the shower front traversing one of the detectors on the ground is 
combined with the times of arrival of light at the photomultipliers of the 
fluorescence telescopes, to enhance significantly the precision with which the 
direction, energy and depth of maximum of the showers can be determined. 
The fluorescence technique has the unique advantage of enabling the energy 
of a shower to be found without resorting to assumptions about hadronic 
physics. This makes it possible to calibrate the shower size found with the 
surface array, and so exploit the near 100% on-time of a surface array. Two 
implementations of the hybrid approach are now operating: the Pierre Auger 
Obse rvatory in Argentina and the Telescope Array in Utah ( Kawai et al.1 
120051 ). Data from the Telescope Array are not yet available and the exposure 
from the Auger Observatory (early 2010) is about 10 times greater. 



4.4 Status of measurements of the energy spectrum 

An energy spectrum based on an exposur e of 12,790 km 2 sr yr has recently 
been r eported by the Auger Collaboration ( Pierre Auger Collaboration et al.1 



l2010bl) . This exposure is a factor 4 higher than achieved with the HiRes de- 
tector at the highest energies and even greater than this at lower energies. 
When deriving the primary energy from fluorescence measurements, an esti- 
mate of the missing energy carried into the ground by hadrons, muons and 
neutrinos must be made based on assumptions about the mass of cosmic rays 
and of the hadronic model. For a primary beam that is a 50/50 mixture of 
protons and iron, it has been found from simulations of showers with the 
QGSJET01 model of hadronic interactions that the energy found fro m the 
fluorescence signal should be increased by 10% ( Barbosa et al.ll2004f h The 



systematic uncertainties on the energy scale sum to 22% with the largest 
coming from the absolute fluorescence yield (14%), the absolute calibration 
of the fluorescence telescopes (9%) and that due to the reconstruction method 
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Fig. 7 The Auger energy spectrum. The spectrum is fitted with two functions [see 
IPierre Auger C ollaborat ion et aLl d2010bf> for details] and the data are compared 
with the stereo spectrum of HiRes (|Abbasi et al.|[2009l ). The uncertainty of the flux 
scaled by E 3 arising from the systematic uncertainty in the energy scale of 22% is 
indicated. 



of the longitudinal shower profile (10%). Efforts are underway to reduce these 
uncertainties. 

The Auger spectrum reported in lPierre Auger Collaboration et all ( 2010fi ). 
and based on over 36,000 events, is shown in Fig. [T] together with the final 
measurement from the HiRes instrument ( Abbasi et al.ll2009t) . It is custom- 
ary to represent parts of the spectrum with a differential form, J = kE 1 . 
It is seen that there is clear evidence of a flattening of the slope at an en- 
ergy of log.E(eV) = (18.61 ± 0.01) where 7 changes from (3.26 ± 0.04) to 
(2.59 ±0.02). At log£(eV) = (19.46 ±0.03) the slope increases very sharply 
to (4.3±0.2). This suppression of the flux is what would be expected because 
of the GZK-effect and is significant at the level of 20 standard deviations. 

An intriguing question is whether the slopes in the region of the suppres- 
sion are different in the Northern and Southern Hemispheres. This would be 
expected if the source distribution was different in the two hemispheres with 
the region with sources at a closer average distance being expected to show 
a less suppressed slope, assuming that the mean injection spectra are the 
same. However the present measurements of (5.5±1.8, HiRes) and (4.3±0.2, 
Auger) arc not significantly different and as the HiRes instrument is not 
longer operating the possibility of any difference will need to be explored 
with future instruments. 

In any case, one may expect that the flux at the very highest energies, 
namely in the GZK cutoff region itself, will be dominated by only a few 
sources within the horizon. For this reason, the assumption of a uniform 
source distribution usually made to produce the synthetic spectra used to 
fit the observed spectrum should not be valid anymore. Even with infinite 
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statistics, the overall spectrum around 10 20 eV will not be of much interest, 
as it will be shaped mostly by the most luminous sources that just happen 
to be here or there in the local sky - in other words, cosmic variance is the 
dominant factor in the region of the spectrum. 

The flux above 10 20 eV, as estimated from the 3 Auger events and 1 HiRcs 
event, is (2A±\%) x l(T 4 km" 2 sr 1 yr 1 or about 1 per square kilometre per 
millennium. It is really quite remarkable that the energy of such rare events 
can be measured to within about 20% and indeed estimates of the energies 
of particle initiating showers at the highest energies are more certain than at 
10 15 eV where there has to be greater reliance on models of hadronic physics. 



4.5 Current status of measurements of the arrival direction distribution 

From a quick look at the first data recorded with the Auger Observatory at 
the highest energies it was evident that no strong sources stood out. It is not 
clear what classes of objects are likely to provide the necessary acceleration 
sites so a well-defined search strategy to look for signals was adopted. Ac- 
tive galactic nuclei, radio-galaxy lobes and gamma-ray bursts have all been 
suggested as possibilities and if these are anisotropically distributed and are 
not too numerous within the distance associated with the GZK-supprcssion 
then an anisotropy might be expected. However there is no prediction for 
source features as compelling as was that of the GZK-cffect for the spectral 
shape. Faced with this dilemma the Auger Collaboration made the decision to 
search for correlations on different angular scales with objects in a particular 
catalogue using energy and distance as addition al parameters. The catalogue 
chosen was the IVeron-Cettv and Veronl ( 20061 ) catalogue (VCV) of quasars 



and active galactic nuclei. Using data recorded between 1 January 2004 and 
26 May 2006 a scan was made for a minimum in the probability P for a 
set of N events expected from an isotropic flux to contain k or more events 
at a maximum angular distance SP from any object in the VCV catalogue. 
P is given by the cumulative binomial distribution X^Lfc C^pi{l — p) N ~ J , 
where p is the fraction of the sky (weighted by exposure) defined by the re- 
gion at angular separation less than tp from a selected source. The parameter 
space was defined by the angular separation ip, the maximum red-shift z max 
and the threshold energy E t h- A minimum in P was found with ip = 3.1°, 
z max = 0.018 (~ 120 Mpc) and E th = 5.6 x 10 19 cV. With these values, 
12 events among 15 correlated with AGNs in the catalogue, whereas only 
3.2 would have been expected by chance. These numbers were found after a 
search and so no probability could be attached to the significance of the asso- 
ciation. However the observation led to the definition of a test to validate the 
result with an indepen dent data set us i ng these pa rameters a priori. Details 
of the test are given in lAbraham et al 1 (120071 120081 ). The test was applied to 
data collected between 27 May 2006 and 31 August 2007 using exactly the 
same reconstruction algorithms as for the initial period. In this independent 
data set, which had almost exactly the same exposure, 8 of 13 events above 
the same energy threshold were found to correlate with events in the VCV 
catalogue. The probability of this occurring by chance is 1.7 x 10~ 3 and the 
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Fig. 8 The directions, in galactic coordinates, of 69 events above 5 x 10 eV 
recor ded at the Pierre Auger Observatory up to 31 December 2009 (|Abreu et alj 
l201Ch . The events are shown as black dots on an Aitoff-Hammer projection with 
the galactic plane across the centre of the diagram. The shaded circles of 3.1° are 
centered at the positions of 318 AGNs in the VCV catalogue that lie within 75 Mpc 
and are within the field of view of the Observatory. The darker shading corresponds 
to greater relative exposure: the exposure weighted-fraction covers 21% of the sky 



result was reported (| Abraham et al.ll2007j) . This correlation corresponds to 
an association of (69j^g)%. 

The Auger Collaboration has continued to take data and as at 31 De- 
cember 2 009 a total of 69 events above the threshold energy have been 
recorded ( Abreu et alJl20ldh . The correlation with objects in the VCV cat- 
alogue, excluding the 15 events used to define the correlation parameters, 
is now (38 ± 6)% with 21 events correlated. The fraction expected from an 
isotropic distribution of event directions is 0.21. From an isotropic flux the 
cumulative probability of such a correlation is P = 0.003. A sky map using 
an Aitoff-Hammer projection is shown in Fig. [5] With the current estimate 
of the correlation at 38%, a 5 sigma signal (P < 3 x 10~ 5 ) will require 110 
events which should be obtained with about 2 further years of operation. 

Using a cross-correlation study the Auger Collaboration have a lso noted 

that t here is an excess of event s from the region around Centaurus A ( Abraham et al.1 
l2007l . l2008llAbreu et al]|2010l) with 13 of the arrival directions within 18° of 
Centaurus A forming 6 pairs separated by less than 4° and 28 pairs separated 
by 11°. Clearly this is a region of sky that will be studied carefully but as the 
search was made a posteriori no statistical significance can be attached to the 
observation. No similar object lies as close to us in the Northern Hemisphere. 

The HiRes collaboration does not find a similarly correlat ion of their 
events with AGNs in the VCV catalogue (jAbbasi et al.ll2008bl) . They find 
that 2 events out of 13 are correlated with VCV objects: this compares with 
5 that would be expected assuming that the correlation is really 38%. The 
difference between 2 and 5 does not rule out a 38% correlation in the northern 
hemisphere, the part of the sky dominantly observed with the HiRes detector. 
It may be that the source distributions arc different in the two parts of the 
sky and, with the steeply-falling spectrum, there may be a small difference 
in the definition of the energy threshold. The situation is inconclusive and 
will not be clarified until new Observatories are built. 
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However that may be, keeping in mind that the excess of correlation 
observed by Auger does not necessarily imply that the AGNs of the catalog 
are related with the sources, and may simply reflect the accidental correlation 
between these AGNs and the actual sources in that part of the sky (convolved 
with the deflections in the intervening magnetic fields), one must be prepared 
to a situation where other sources in the Northern hemisphere lead to a 
different level of correlation with the AGNs of that hemisphere. In other 
words, there is no direct conflict between the observations currently available, 
and a systematic study of the entire sky with comparable data sets should 
tell us more about the origin of the detected anisotropy, and what it implies 
concerning the UHECR sources and deflections. 

In the mean time, it is interesting to note that the current data sets do not 
reveal obvious multiplets or close correlations with known sources. This may 
be due to a very large number of sources, from which only one event at most 
has typically been observed so far, or to relatively large deflections, causing 
the sources to overlap in the sky. In both cases, the observational solution 
consists in increasing the statistics at high energy, where the deflections are 
reduced and the horizon scale becomes smaller, so that only a few sources can 
significantly contribute. The fact that some anisotropy is already detectable 
gives strong confidence that the expected breakthrough is within reach of the 
next generation detectors, pro vided they can gat her a statistics of the order 
of a few 10 5 km 2 sryr or more (|01into et aLll2009f) . 

From the moment when the first few significant hot spots are observed in 
the UHECR sky, a new era will begin, with the study of individual sources. 
As indicated above, this should lead to key parameters such as the source 
power and spectral shape. The current statistics does not allow one to confirm 
whether the excess of events possibly seen around the location of Centaurus 
A is related to the first of such sources, but if it is the case, then the next 
generation of detectors should accumulate tens or even hundreds of events 
from that or similar sources. 



4.6 Measurement of the chemical composition of the primaries 

To interpret the information that we have on the energy spectrum and arrival 
direction pattern of ultra high-energy cosmic rays we need to know the charge, 
i?, of the particles, ideally on an event- by- event basis. Then, for example, if 
one had identified a point source, the bending of the particles from it as 
a function of energy could be used to distinguish between different models 
of the magnetic fields lying between the source and earth. However it is 
very difficult to be sure even of the atomic mass, A, as the only methods 
that are practical for estimating it require assumptions about the hadronic 
interactions at centre-of-mass energies well above what will be reached with 
the LHC. 

The concept of the most promising method is to compare the average 
depth of shower maximum (called A max ) as a function of energy for pro- 
tons, iron nuclei and photons with what is predicted for different hadronic 
models. The depth of maximum can be determined to ~ (20 30) g/cm 2 
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Fig. 9 Left: The average measurements of X max as a function of energy reported 
by the Auger Collaboration b ased on 3754 events c ompared with the predictions 
of 4 different calculations (sec Abraham et al. 2010|, for details). Right: The mea- 
surements of the rms of -X max as a fu nction of energy for t he data of the figure on 
the left, compared with calculations (jAbraham et al.H2010l ). 



using the fluorescence technique. For a mass composition of a single species, 
I max is expected to increase with energy as 7-rays from tt° decay become 
more energetic and so the position of shower maximum is pushed deeper into 
the atmosphere. The average A" max for iron nuclei lies below that for pro- 
tons because the energy per nucleon of an incoming primary is smaller for 
heavier nuclei at a given energy. The change of X max with energy is known 
as the elongation rate. In addition to the elongation rate being a useful in- 
dicator of primary mass, a key parameter is the fluctuation in the depth of 
shower maximum at a particular energy. It is intuitively obvious that as an 
iron nucleus is larger than a proton, the cross-section for interaction will be 
greater. Thus fluctuations in the position of shower maximum also contain 
information about the primary mass. However, although these parameters 
can be measured rather accurately, the interpretation of the measurements 
must rely on models of shower development. 

Wo rk on this problem ha s recently been reported by the A uger Collabo- 
ratio n ( Abraham et al.|[2010[ ) and by the HiRes Collaboration ( Abbasi et al.1 
l2010h . With the Auger Observatory the depth of maximum has been de- 
termined for over 3700 events of energy > 10 18 eV. Each event has been 
reconstructed using the hybrid method: 2% of them are stereo events for 
which two independent estimates of X max are made. Using stereo informa- 
tion it can be demonstrated that the accuracy of estimation of X max in a 
non-stereo event is ~ 20g/cm -2 (about 200 m at the height of the maximum 
of a typical shower, 750g/cm~ , or 2.5 km above sea-level in the vertical 
direction). 

The Auger data on the depth and width of the fluctuation of shower 
maximum are shown in Fig. [5J The predictions for proton and iron nuclei 
of AT max and rms(J max ) are given for different models of the hadronic inter- 
actions (QGSJET01 etc). It is clear that if these models are approximately 
correct then the mean mass of the cosmic rays is becoming heavier as the 
energy increases. However the results from the HiRes collaboration seem to 
tell a different story. They conclude that their data, totaling 815 events, is 
entirely consistent with proton primaries from ~ 3 x 10 18 eV. The reasons 
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for the differences between the two measurements are not understood. The 
HiRcs result is highly dependent on Monte Carlo calculations whereas the 
Auger measurement makes very little use of such calculations except, as is 
true for both measurements, when a comparison is made with the predic- 
tions. The Auger Collaboration have taken precautions in their analysis to 
test that the effects observed are not dependent on zenith angle, demon- 
strating that deeply penetrating events are not missed, and it is hard to see 
how distributions can be made too narrow as a consequence of the analysis 
procedures. 



5 Some conclusions and perspectives 

Magnetic hclds play a key role in the acceleration, propagation and radiation 
of galactic cosmic rays. Current generation of X-ray Chandra, XMM-Newton 
and Suzaku and gamma ray telescopes H.E.S.S., MAGIC, MILAGRO, VER- 
ITAS and Fermi have provided conclusive evidences for CR particle acceler- 
ation up to the energies at least about 100 TeV in young supernova remnants, 
starburst regions, radio galaxies and AGNs . With the coming gen eration of 
the fo cusing hard X-ray tel escopes Astro-H ( Takahashi et al.ll2010l) and NuS- 



TAR (jHarrison et alJ2010h new possibilities will be opened to study in details 



hard X-rays images and spectra that of cosmic ray sources providing unique 
information on both particle spectra and fluctuating magnetic fields in the 
sources. A project of very large X-ray telescope IXO that is currently under 
discussion in addition to the un precedent sensitivity will be likely supplied 
with X-ray imaging polarimetry (jBarcons et al.ll201l[ ). Imaging polarimetry 
is sensitive to details of magnetic field amplification mechanism i n the particle 



accel eration sources like supernova remnants and radio galaxies (jBvkov et al.1 
l2009h . The Cherenkov Telescope Array CTA project under development is 
the next generation ground-based gamma-ray instrument, which is supposed 
to provide at least 5 times better sensitivity in the current energy domain of 
about 100 GeV to about 10 TeV and an extension of the accessible energy 
range well below 100 GeV and to above 100 TeV. 

At extremely high energies, E ~ 10 20 eV, the impact of galactic and ex- 
tragalactic magnetic fields on the propagation of cosmic rays becomes less 
dramatic, which should result in large and small scale anisotropics of the 
cosmic rays. Although the corresponding fluxes are very low, considerable 
progress has been made in the recent years, with the detection of the giant 
cascades (so-called extensive air showers) induced by these cosmic rays in 
the atmosphere, by large detectors consisting of ground arrays (AGASA), 
fluorescence telescopes (HiRcs) or hybrid detectors (Auger, currently in op- 
eration) . 

The extension of cosmic ray studies up to energies of the order of 10 20 eV 
is likely to give rise to what is referred to as "proton astronomy" (or more 
generally "charged particle astronomy" ) , where charged cosmic rays are used 
as astronomical messengers pointing roughly back to their sources, because 
of the reduced deflections in the intervening magnetic fields. Another key 
aspect of the cosmic rays phenomenology in this energy range is the exis- 
tence of a physical horizon, associated with the energy losses suffered by the 
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cosmic rays as they interact with the background radiation fields (notably 
the cosmological microwave background). 

Whether protons or heavier nuclei, particles of such high energies can 
arrive only from relatively nearby accelerators, say within 100 Mpc, which 
decreases dramatically the number of potential sources capable of acceler- 
ating cosmic rays up to > 10 20 eV and contributing to the observed flux. 
In principle, many weak sources may be at work, which may prevent one 
from detecting them individually, by the accumulation of events in a given 
direction. However, even with large source densities, a few most nearby and 
brightest sources may be expected to contribute a significant fraction of the 
total flux, in a similar way as what is observed at all photon wavelengths, 
where a few sources dominate the sky. Alternatively, there may be only a few 
sources within the horizon, especially if one considers the tough req uirements 
associated with 10 20 eV proton accelerators ( Aharonian et all 2002h . This ex- 



cludes, in particular, astrophysical objects like ordinary galaxies, unless these 
galaxies accelerate the highest energy cosmic rays through transient events 
related to compact objects like Gamma Ray Bursts. Whether we can iden- 
tify the accelerators of extragalactic cosmic rays using the highest energy 
protons is a question which largely depends on the strength and structure 
of the large scale IGMF. However, if the latter are not too large, as seems 
to indicate the current UHECR data showing anisotropy above 55 EeV or 
so, the cosmic rays from identified sources may be used in turn to probe the 
intergalactic magnetic fields in the nearby universe, as well as the structure 
and strength of the galactic magnetic field, from the observed (energy- and 
charge-dependent) deflections. 

The identification of the first cosmic ray sources by astronomical means, 
i.e. through hot spots or elongated multiplets with a characteristic rigidity 

ordering, may be provided by the next generation detectors, either on the 

groun d, with the Northern site of the Pierre Auger Observatory (jBlumer and the Pierre Auger Collaboration] 

l2010f ). proposed for construction over 20,000 km 2 , with the proven technique 

of Auger South and similar pe rformance, or in space, with the JEM-EUSO 

mission ( Ebisuzaki et al.l I2010T ) . proposed to be installed on the Japanese 

module of the International Space Station, for an even larger aperture at the 

highest energies. This is expected to lead to a new era of cosmic ray studies, 

where nearby sources may be studied individually, and particle acceleration 

may be explored in the most challenging regime, up to above 10 20 eV. 

Finally, we stress that the strong coupling between cosmic ray studies 
and the understanding of the magnetic field generation, amplification and 
structure in the universe calls for an increased exchange of knowledge between 
the different communities involved. Important progress is expected in both 
fields in the coming years, and any progress in one of the various aspects 
of the problem will benefit to the others, and to the whole field of high- 
energy astrophysics, where energetic particles are responsible for non-thermal 
emission at all energies, as well as to the physics of the interstellar medium 
in general, in which the magnetic field is a key component. 
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